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Abstract. Monte-Carlo implementations of /cr-factorization formula with both KLN and 
q | running-coupling BK unintegrated gluon distributions for nucleus-nucleus collisions are used to 

, analyze recent experimental data on the particle multiplicities from RHIC(Au+Au@200GeV) 

£^ ' and LHC(Pb+Pb@2.76TeV). We also compare the predicted transverse energy at midrapidity 

to new data from ALICE. 
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Q\ . 1. Introduction 

Relativistic Heavy Ion Collider (RHIC) at BNL has started in 2000 and many interesting data 
have been obtained. Experiments at Large Hadron Collider (LHC) at CERN has just started. It 
is well known fact that nearly perfect fluid picture works to explain large elliptic flow discovered 
at RHIC and LHC. However, there is no detailed understanding on the particle production in 
high energy hadronic collisions. Especially subsequent non-equilibrium evolution of the system 
(Glasma) created in nuclear collisions is not known well. Without knowing the detailed glue 
dynamics in early stages of nucleus collision, correct initial conditions for hydrodynamics would 
not be obtained. First important issue is to understand the gluon production in the first moment 
of collision. In order to study such gluon production in nucleus-nucleus collision, we shall use 
the Monte-Carlo implementation of /^-factorization formulation in which fluctuations of the 
position of nucleon inside a nucleus is taken into account, and we can study nucleus-nucleus 
collision event-by-event. 

2. Theoretical Models 

We will use the Monte-Carlo KLN (MC-KLN) model [1] and its extension to running coupling 
Balitsky-Kovchegov (MCrcBK) [2] for the computation of gluon production in heavy ion 
collisions. Both model apply the &t-factorized formula [3] in the transverse plane perpendicular 
to the beam axis locally. The number distribution of produced gluons is given by 

dN g AN r 



/d p_\_ f 
J d 2 k± a s (t> A (xi,k 2 ± ) (pB(x2,(P±- k ±) 2 ) > (!) 
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with N c = 3 the number of colors. Here, p± and y denote the transverse momentum and the 
rapidity of the produced gluons, respectively. The light-cone momentum fractions of the colliding 
gluon ladders are then given by x± t 2 = p± exp(±y)/ \/siyN, where ^/snn denotes the center of 
mass energy. A constant gluon multiplication factor n g = 5 is assumed in the MCrcBK model 
to obtain pj_-integrated hadron yields [2] . It does not appear in the corresponding formula for 
the transverse energy dE±/dy since gluon splitting and hadronization conserves energy. 

At each grid point, we compute the thickness function Ta(t±) to obtain the local saturation 
scale and then compute gluon production probability. Within a hard disk nucleon approximation, 
thinkness function is given by 

„, , , number of nucleon within S , . 

T A (r±) = ^ (2) 

where we assume that the area S has the same value as the inelastic proton-proton cross section 
at the incident energy of y / i]v7V = 200 GeV independent colliding energy. 
In MC-KLN model, saturation momentum is defined as 

Q ^ ;r± ) = 2G e V*-^^(^)\ (3) 
1.53 fm V x J 

where A is a free parameter which is expected to have the range of 0.2 < A < 0.3 from HERA 
global analysis. In MC-KLN, we assume the gluon distribution function as 

^^"Jw—mU)' <4) 

On the other hand, in MCrcBK, is obtained from the Fourier transform of the numerical 
results of the running coupling BK (rcBK) evolution equation [7J[8l[9]: 

4>A,B{k,x,b) = „ C f n .„ / d\ e- ik - r V r 2 AA G (r,y = ln(x /x),6). (5) 



a s (k) (2tt) 3 

where Cf = (N% — 1)/2N C , xq = 0.01, and Mg is related to the quark dipole scattering amplitude 
that solves the rcBK equation, N as follows: 

Af G (r,x) =2Af(r,x) - Af 2 (r,x) . (6) 

In MC-KLN model, x dependence is determined by eq. ([3]), but in rcBK, x dependence can 
be obtained from the equation. Therefore, we expect that MCrcBK model has more predictive 
power than MC-KLN model. 



3. Results 

3.1. MC-KLN 

In the figure [H the charged particle multiplicities at RHIC and LHC as a function of -/V part from 
KLN, fKLN, and MC-KLN model are compared to PHOBOS data. 
In the KLN model, we plot the equation 



1 dN 
A/part dy 




(7) 



This equation is obtained by picking up the most important contribution of integration in kt- 
factorized formula by using average nuclear saturation momentum. On the other hand, fKLN 
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Figure 1. Centrality dependence of charged 
hadron multiplicity in Au+Au collision from 
the KLN, fKLN, and MC-KLN model are 
compared with PHOBOS data [5] 
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Figure 2. Centrality dependence of charged 
hadron multiplicity in Pb+Pb collision at 
^/i^T = 2.76 TeV from the fKLN, and MC- 
KLN model are compared with ALICE data 
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model [1] does not take into account the effects of fluctuations of nucleons inside a nucleus unlike 
the MC-KLN model. The main difference between KLN and fKLN model is that fKLN model 
uses the local saturation scale instead of assuming a average nuclear saturation momentum as in 
KLN. Therefore, we see that computations with local saturation momentum improve results. A 
comparison between fKLN and MC-KLN model shows the effects of fluctuations of the position 
of nucleon inside a nucleus. One sees such effect in the peripheral collision which suggests that it 
is important to take into account fluctuation effect in the discussion of the centrality dependence 
of particle multiplicities. 

In Fig. [21 the centrality dependence of charged hadron multiplicity in Pb+Pb collision at 
y/ s NN = 2.76 TeV from ALICE experiment is compared to the results from fKLN and MC- 
KLN. Where we use A = 0.23 in Eq. ([3]) which controls the x-evolution speed. If A is small, 
multiplicity becomes small, and we may think additional particle production mechanism during 
the evolution of the system, especially before thermalization. 



3.2. MCrcBK 

We now compare in Figj3l the centrality dependence of charged hadron with the results obtained 
from the MCrcBK model in which unintegrated gluon function is taken from the numerical 
solution of the rcBK equation. The initial condition for the rcBK evolution is assumed to be 
MV model: 



A/"(r,y = 0) = 1 - exp 
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with the initial saturation scale Qso 
coupling, we use 



= 0.2 GeV 2 for nucleon and A 

f 2, 127T 

a s [r ) 
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4C 2 
r 2 A 2 



(8) 

0.2 GeV. For the running 
(9) 



with C = 1 in solving rcBK equation. rcBK unintegrated gluon function with these parameter 
set describes the centrality dependence of charged hadron multiplicity for both RHIC and LHC 
which indicates that most of the gluon is produced at the first impact of the nuclear collision. 
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Figure 3. Centrality dependence of charged hadron multiplicity in Pb+Pb collision at 
^snn = 2.76 TeV and Au+Au collision at ^s^n = 200 GeV from the MCrcBK model are 
compared. 



However, we should check parameter dependence carefully. For example, if one changes the 
parameter C in the running coupling, evolution speed changes. 

In the bottom up scenario [TO], the number of gluon produced just after the collision will 
increase by a factor of a s (Q 2 )~ 2 / 5 during the subsequent evolution of the system. It would be 
interesting to include such effect into MCrcBK calculations. 

3.3. Gaussian shape 

So far we assumed that nucleus-nucleus collision is described by the incoherent sum of nucleon- 
nucleon collisions which will occur when transverse distance squared between two nucleons is 
smaller than the inelastic proton-proton cross section &nn divided by it: 

0*-*j) a + (w-w) 2 <^ (10) 

This amounts to assume that nucleon is hard sphere (or disk). However, this approximation 
may not work in very high energy hadronic collisions. Let us think about the Gaussian shape 
nucleon in order to take into account the effect of extension of nucleon size as incident energy 
increases. In this case the thickness function becomes 

Tp(r)= 2^ eXp[ " r2/(2jB)] ' 
and the probability of nucleon-nucleon collision P(b) at impact parameter b is 

P(b) = 1 - exp[-kT pp (b)}, T pp (b) = J d 2 sT p (s)T p (s - b) . (12) 

where (perturbatively) k corresponds to the product of gluon-gluon cross section and gluon 
density squared. We fix k so that integral over impact parameter becomes the nucleon-nucleon 
inelastic cross section a^N at the given energy: 



a NN = / d 2 b(l - exp[-kTpp(b)]) 



(13) 



In this work, we use B = 0.2 fm 2 , ctnn = 41.94 mb for y / J/v r /v — 200 GeV and ct^n = 61.36 mb 
for -y/F/vTv = 2.76 TeV. The result of this model is plotted in right hand side of Fig. [3j Model 
underpredicts the multiplicity at peripheral collisions. One possible interpretation may be the 
following: viscosity is large at large impact parameter, entropy production may be larger as 
impact parameter becomes large. 



3-4- Eccentricity and transverse energy 
Finally, we plot eccentricity defined by 



(y 2 



(y 2 



(14) 



in Fig. HI Since the eccentricity is proportional to the magnitude of elliptic flow, it is important 
to know the initial value of this quantity. One sees that the incident energy dependence of the 
eccentricity is very small according to the results from MCrcBK. 
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Figure 4. Centrality dependence 
of eccentricity with respect to 
reaction plane in Pb+Pb collisions 
at y/sjvN = 2.76 TeV and in 
Au+Au collisions at ^/snn = 200 
GeV. 




Figure 5. Centrality dependence of the transverse energy in Au+Au collision at tJsnn = 200 
GeV (left) and Pb+Pb collision at ^/i^F = 2.76 TeV (right). PHENIX data from PH, ALICE 
data from [T2] , 



Fig. [5] shows the centrality dependence of the transverse energy at central rapidity for the 
MV initial condition, and for an initial condition featuring a more rapid fall-off at k± > Q s ; 
see ref. [2] for details. We note that ~ 2.5% (0.5%) of the energy of the beams is predicted 
to be deposited initially into the central rapidity region in central collisions. Longitudinal 
hydrodynamic expansion (— pAV work) may reduce the transverse energy by up to a factor 
of 2 |13j . In all, the model based on rcBK evolution of the gluon distribution appears to be 
consistent with the multiplicity and transverse energy data over a good order of magnitude in 
both N P art and y/s NN . 



4. Summary 

Within the MC-KLN and MCrcBK model, we compute the gluon production based on the 
picture of Color Glass Condensate. Both MC-KLN and MCrcBK model reproduce the centrality 
dependence of charged hadron multiplicity at RHIC and LHC energies which indicate that the 
entropy production during the thermalization process may be small as well as the small viscosity 
after thermalization. However, we need detail systematic study of particle production by looking 
at different observables together with the parameter dependence. 
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